function [y] = fun_weifen1(a,b)
y=-b*(1+a*b)%%y=-b*(1+a*b)  
end

